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ABSTRACT 

The 11 Mpc Ha and Ultraviolet Galaxy (11HUGS) Survey traces the star formation 
activity of nearby galaxies. In addition within this volume the detection completeness 
of core-collapse supernovae (CCSNe) is high therefore by comparing these observed 
stellar births and deaths we can make a sensitive test of our understanding of how stars 
live and die. In this paper, we use the results of the Binary Population and Spectral 
Synthesis (BPASS) code to simulate the 11HUGS galaxies’ Ha and far-ultraviolet 
(FUV) star formation rate indicators (SFRIs) and simultaneously match the core¬ 
collapse supernova (CCSN) rate. We find that stellar population including interacting 
binary stars makes little difference to the total CCSN rate but increases the Ha and 
FUV fluxes for a constant number of stars being formed. In addition they significantly 
increase the predicted rate of type Ibc supernovae (SNe) relative to type II SNe to 
the level observed in the 11HUGS galaxies. We also find that instead of assuming a 
constant star formation history (SFH) for the galaxies our best fitting models have a 
star formation rate (SFR) that peaked more than 3 Myrs ago. 
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1 INTRODUCTION 


Core-collapse supernovae (CCSNe) are the explosive 
events t hat mark the death of massive stars (> 8 Mq) 
dSmartt et alj[2009h . They are classified into hydrogen-rich 
events, Type IIP/L/n/b, and hydrogen-deficient events, 
type Ib/c, according to th e presence or ab sence of hydrogen 
lines in their spectrum (IFilippcnkol 1199711 . This is directly 
related to the composition of the outer layers of a progenitor 
star at the moment of explosion. Type II supernova (SN) 
progenitors have kept their hydrogen rich outer envelope 
intact, while type Ib/c have lost the hydrogen envelope. 
The most direct evidence linking the star that explodes to 
the type of SN is from direc t detection of the progenito r 
in pre-explosion images (e.g. iLi et al . ISqM mar ia. 
Non-detections of progen itors i n pre-exp losion images for 
type Ibc SN were used by [Eldridge et al.1 (|2013|) to indicate 
the important effect of interacting binary stars on determin¬ 
ing the CCSN progenitor diversities and local CCSN rate. 
This idea was later reinforced by the probable detection of 
a binary system progenitor for the type Ibc SN iPTF13bv n 
llCao et al.ll2013l : lBersten et al.ll2014 : [Eldridge et ahll2015l l. 


However, it is difficult to directly detect CCSN pro¬ 
genitors because the event must occur in a galaxy within 
about 30 Mpc that has been imaged by the Hubble 
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Space Telescope before explosion. Therefore indirect 
methods must be exploited to gain further insight into 
CCSNe. The most useful methods include the relative 
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One other interesting constraint on CCSNe and their 
progenitors can be obtained by exploiting the link between 
stellar deaths a n d ste llar births. This was attempted by 
iBotticella et al.l (I2012T ). They used the 11HUGS galaxy 
catalogue that lists the Ha and FUV fluxes for a large 
number of galaxies within 11 Mpc. These fluxe s are com- 
monly used tracers of recent star formation dKennicuttl 
Il998l ). The amount of Ha and FUV emitted are determined 
by the most massive hot stars that evolve quickly and 
exp lod e as CCS Ne a fe w million years after their formation. 
IBotticella et al.l J2012ih were able to show that the observed 
Ha and FUV fluxes from the 11 HUGS galaxies were 
consistent with the observed CCSN rate for the galaxies if 
the minimum stellar mass for a CCSN to occur was 8Mq, 
in a greement wi th res ults from CCSN progenitor detection 
in ISmartt et al.l (l2009h . 
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This study was repeated bv lHoriuchi et all J 2 OI 3 I ) who tried 
a slightly more sophisticated modelling of the 11HUGS 
galaxies. Rather than using the canonical expected SFR 
conversion factors from iKennicuttl (1 19981) they estimated 
the expected Ha, FUV flux and CCSN rate from the 
l atest stel lar models including rotation calculated by 
lEkstrom et al.l d2012l '). Over the last few decades the 
understanding of how rotation affects st ellar evolution has 
signif i cantly advan ced (see reviews by iMaeder fc Mevnetl 
120001 : iLangerd 120121 ') . The general effect is to increase the 
main-sequence life time and temperature of massive stars 
as they evolve so lead to m ore Ha a nd FUV flux being 
emitted for the same SFR. iHoriuchi et aid d2013l 'l found 
that it is possible that there could be a CCSN rate excess 
as well as that the observations disfavoured the rotating 
models as they made this SN rate excess worse. 


There is a second important factor left for our stellar 
population models that is the effect of interacting binaries. 
When two stars are orbiting each other in a binary-star 
system the opportunity for mass transfer between the 
two stars arises. The more massive star typically evolves 
faster than the less massive star, fills its Roche-Lobe and 
transfers mass to the companion star. This produces more 
hot and luminous stars and therefore has simil ar observed 
consequence as stellar rotation does (e. g. Vanbeveren et all 
1 19981 : lEldridge et al.l [20081 : IZhang fc Lill2012h . However the 
situation may be more complex as there is evidence that 
the most rap i d rot ation only occurs in binary systems 
(Ide Mink et al.l 12013 :1. 


In th is pa per we continue the work of iBotticella et al.l 
ll2012l ) and IHoriuchi et, ~ahl (|2013l ) by attempting to match 
the observed Ha, FUV fluxes and CCSN rate in the 11HUGS 
galaxies. We will use the Binary Population and Spectral 
Synthesis (BPASS, http://bpass.auckland.ac.nz) code 
that accounts for the impact of interacting binaries on 
stellar populations (Eldridge & Stanway 2009, 2012). The 
BPASS code has been tested against a number of different 
observables that are dependant on the evolution of massive 
stars. Our aim is to apply a simple test to BPASS models to 
reproducing the CCSN rate, Ha and FUV fluxes. Our work 
concentrates on the effect of interacting binary stars on the 
CCSN rate of their host galaxies. We also include two other 
refinements to our models. First we reproduce the chemical 
compositions of the 11HUGS galaxies based on the B-band 
magnitude of nearby galaxies sample rather than using a 
Solar composition population. Second we allow the SFH to 
vary rather than assuming only a constant SFR. 


The structure of this paper is arranged as follows. In 
Section 2, we introduce the 11HUGS and CCSN samples 
and how we determine the metallicity of the 11HUGS 
galaxies. In the following Section 3, we summarise the effect 
of interacting binary stars on CCSN rate and SFRIs we 
find from the BPASS code’s results. We also outline how 
we use these BPASS results to construct our synthetic 
galaxy models. In Section 4, we compare our models 
to observations before finally discussing our results and 
presenting our conclusions in Section 5. 


Table 1. Nearby g alaxy sample from t he 11 HUGS sur¬ 
vey llLee et_al]j20n^_ [Botticella et ahl 120121) and CCSNe sam¬ 
ple from lEldridge"et al. | d2013i) : Numbers of galaxies, the to¬ 
tal Ln a /10 43 ergs —1 , Luv/10 28 ergs — 4 Hz —1 and CCSN rate 
(Rcc/yr -1 )- 1 The T.i y for Sample A is an estimation with re¬ 
spect to the relative ratio of Ha and UV flux from Sample B and 
C. 


Parameter 

Sample A 

Sample B 

Sample C 

Ng a l 

383 

312 

167 

L Ha /10 43 ergs — 1 

1.1 ±0.05 

0.9 ± 0.05 

0.7 ±0.04 

L 0 v/lO 28 ergs _ 1 Hz _1 

Moots 

88 ±6 

67 ±4 

Ncc 

17 

15 

14 

NsnI 

4 

4 

3 

NsnII 

13 

11 

11 

Rcc/yr -1 

1.13 ±0.27 

1 ± 0.26 

0.93 ±0.25 


2 GALAXY AND CCSNE SAMPLES 

The galaxy samples w e adopt in this paper are drawn from 
IBotticella et al.l (I2OI2I) . There are three different samples: 
383 galaxies with measured flux in Ha (sample A), 312 
galaxies with measured fluxes in Ha and FUV (sample B) 
and 167 galaxies with measured fluxes in Ha, FUV and 
Total Infrared (TIR) (sample C). We list the observed 
fluxes of these samples in Table |T] We broadly assume that 
the relative ratio of Ha and FUV flux is constant. Therefore 
applying the averaged flux ratio of Ha and FUV of Sample 
B and C to Sample A, we estimate the FUV luminosity of 
Sample A: Luv = (106 ± 8) x 10 28 ergs -1 Hz -1 . 

In addition, we list the details of the 17 CCSNe dis¬ 
covered within 11 Mpc from 1998 to 2012 in Table [2] It is 
directly based on the work of IBotticella et ahl ((20121 ) . where 
they picked out all identified SNe occurred in galaxy sample 
A of 11HUGS s urvey from 18 85 to 2010 from the Asiago 
SNe catalogue dBarbon et al] 120081) . We limit our SNe 
sample to only core-collapse events. We al so exten d this 
CCSN list up to 2012 using the SN list from lEldridge et al.l 
(|2013h . This list of SNe was compiled by searching the 
catalogue maintained by the International Astronomical 
Union (LAU) Central Bureau for Astronomical Telegrams 
(CBETjj. It was searched for all SNe discovered in the 
fifteen-year period from 1 Jan 1998 to 30 March 2012. 
The host galaxies of these selected SNe have a recessional 
velocity (corrected for Local Group infall on Virgo) of 2000 
kms -1 or less, corresponding to a distance limit of 28 Mpc 
for Ho = 72kms -1 Mpc -1 . 

All CCSNe used in our work is in the field of 11HUGS 
survey and related to the adopted SFRIs. Of these there are 
13 type II SNe and 4 type Ibc SNe. Among all the CCSNe 
host galaxies, the host galaxy of SN 2005at is in Sample A 
and B but not in Sample C, and the host galaxy of both SN 
2005af and SN 2011ja is only in Sample A. It is important 
to consider whether our sample of SN is complete and 


1 http://www.cbat.eps.harvard.edu/lists/Supernovae.html 
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Table 2. CCSNe discovered within 11 Mpc from 1998 to 2012 
as listed bv lEldridge et all (|2013h . X SN 2005at host galaxy NGC 
6744 is in Sample A and B not in Sample C , 2 SN 2005af and SN 
2011ja host galaxy NGC 4945 is only in Sample A. 


Table 3. Metallicity classification of Sample B galaxies: the frac¬ 
tion of total galaxies, Hq flux, FUV flux and number of CCSNe 
when the 11 HUGS galaxies are divided into the 5 metallicities of 
BPASS. 


SN 

Gal. 

Type 

M b 

D/Mpc 

2002ap 

NGC 628 

Ic 

-19.65 

7.3 

2002bu 

NGC 4242 

Iln 

-18 

7.4 

2002hh 

NGC 6946 

IIP 

-20.65 

5.9 

2003gd 

NGC 628 

IIP 

-19.65 

7.3 

2004am 

NGC 3034 

IIP 

-19.08 

3.53 

2004dj 

NGC 2403 

IIP 

-18.77 

3.22 

2004et 

NGC 6946 

IIP 

-20.65 

5.9 

2005cs 

NGC 5194 

IIP 

-20.7 

8 

2007gr 

NGC 1058 

Ic 

-18.2 

9.2 

2008ax 

NGC 4490 

lib 

-19.39 

8 

2008bk 

NGC 7793 

IIP 

-18.41 

3.91 

2009hd 

NGC 3627 

IIP 

-19.61 

10.05 

2011dh 

NGC 5194 

lib 

-20.7 

8 

2012fh 

NGC 3344 

Ic 

-18.8 

6.6 

1 2005at 

NGC 6744 

Ic 

-20.9 

9.4 

2 2005af 

NGC 4945 

IIP 

-20.6 

4.08 

2 2011ja 

NGC 4945 

IIP 

-20.6 

4.08 


whether there might have been some CCSNe unobserved. 
For example some CCSNe might occur in very dusty regions 
of galaxies and remain unobserved or some SNe might for m 
black holes and remain unobserved (e.g. ISmar tt) 12015 I) . 
However modern SN searches are unlikely to have missed 
many SNe. In addition if any SNe are unobservable due 
to dust, their associated Ha and FUV emission will also 
be decrease thus in this work should not contribute to the 
observed SFRIs. Although in light of these concerns we can 
suggest that our CCSN rate is at least a robust lower limit. 


3 METHOD 

Our method to reproduce the 11HUGS Ha flux, UV flux and 
CCSN rate falls into three steps. First a metallicity classi¬ 
fication of the 11HUGS galaxies is applied so that we can 
blend our models of different metallicities and obtain a more 
realistic simulation of 11HUGS sample. Second we describe 
the BPASS stellar population model that we use to construct 
our models, and discuss the effect of interacting binary stars 
on CCSN rate and SFRIs of star-forming galaxies. Finally 
we outline our method to create a synthetic galaxy model 
with various star-formation histories. 



Z 

[O/H] 

fgal 

fHa 

fFUV 

fcCSN 


0.001 

7.5 

0.178 

0.003 

0.002 

0 


0.004 

8.1 

0.451 

0.038 

0.033 

0 

T04 

0.008 

8.4 

0.197 

0.134 

0.096 

0.057 


0.020 

8.8 

0.117 

0.324 

0.253 

0.531 


0.040 

9.0 

0.057 

0.501 

0.616 

0.412 


0.001 

7.5 

0.083 

0.001 

0.001 

0 


0.004 

8.1 

0.292 

0.009 

0.009 

0 

KK04 

0.008 

8.4 

0.317 

0.052 

0.04 

0 


0.020 

8.8 

0.225 

0.339 

0.258 

0.353 


0.040 

9.0 

0.083 

0.599 

0.692 

0.647 


0.001 

7.5 

0 

0 

0 

0 


0.004 

8.1 

0.498 

0.052 

0.043 

0 

PP04 

0.008 

8.4 

0.419 

0.427 

0.324 

0.294 

03N2 

0.020 

00 

bo 

0.079 

0.521 

0.633 

0.706 


0.040 

9.0 

0.003 

0 

0 

0 


This shows generally good agreement with other local sam¬ 
ples of galaxies. The relation is, 

[0/H\ = 12 +log (O/H) 

= -(0.185 ±0.001)M fl + (5.238 + 0.018) (1) 

The ratio of oxygen and hydrogen [O/H], is a common 
measurement of galaxy metallicity and able to be estimated 
by the B-band magnitude of the host galaxy taken from 
I Lee et all ( 201 1+ However, this L-Z relation is strongly 
affected by dust attenuation and the choice of which 
str ong-line abundan c e cali brations are used as discussed 
in [Kewlev & Ellison] |2008|). To quantify the influence of 
dust, Tremonti et aL (12004 1) correct the galaxy luminosities 
for_ i ntrinsic a t tenua tion using the attenuation curve of 
I Charlot fe Falll (|2000l ) and assume that the stars experience 
one third of the attenuation measured in nebular gas. To 
account for the effect of using different strong-line ab un¬ 
dance calibrations, we use the work of iKewlev fe Ellisonl 
(120081) who derived conversion relations between different 
strong-line calibrations. Here we use their conservation 
relations b etween Tremont i model (T04 ) and two other 
mod els of iKobu lnickv fc Kewlevl (120041 ) model (KK04) 
and IPettini fe Pagell 1 20041) model (PP04 03N2). These 
two models are the possible extremes to the T04 model. 
Therefore by comparing to these we can estimate how our 
results depend on the metallicity classification assumed. 


3.1 Metallicity Classification 

Previous studies have only considered reproducing the stel¬ 
lar populations in 11HUGS samples by using only Solar 
metallicity composition for the stellar populations. However 
it is well known that synthetic population results depend 
on metallicity even when interacting binaries are included. 
Therefore we take account of the varying chemical compo¬ 
sitions of the 11HUGS galaxies. We estimate the galaxies’ 
metallicity by using the luminosity- metallicity (L -Z) r elatio n 
of the SDSS galaxies provided by iTremonti et ali (120041) . 


BPASS has five different metallicities for its model 
populations, ranging from one twentieth to twice Solar. 
To determine how much star formation occurs at different 
metallicity environments we compare the galaxies’ metallic¬ 
ity to the input values for the BPASS models and calculate 
the fraction of total observed galaxy population (f ga i), Ha 
flux (fHa), FUV flux (fFUv) and CCSN rate (fccsN) at 
each metallicity. For this classification we use Sample B 
because it is large enough to cover more than 80 per cent 
of the total sample and has both Ha and FUV observations 
available. Our metallicity classifications with respected to 
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Table 4. Predicted flux of Ha and FUV at SFR of lMgyr 1 , lOMgyr 50Mgyr 1 , lOOMgyr 1 with estimated uncertainty, due to 
stochastic effect of constructing the synthetic galaxy models from individual star clusters whose masses were selected at random. 


SFR 

Z 

IW(10 41 

erg s 1 ) 

Lpuv/(10 28 

erg s 1 Hz 1 ) 

Rccsn/( 10 2 yr) 

Rlbc /Rll 



single 

binary 

single 

binary 

single 

binary 

single 

binary 


0.001 

3.1 ± 1.2 

4.4 ± 1.5 

1.3 ±0.2 

1.8 ±0.2 

1.2 ±0.2 

1.1 ± 0.2 

0.003 ± 0.003 

0.46 ±0.04 

IMg/yr 

0.004 

2.5 ±0.9 

5.7 ± 1.3 

1.1 ±0.2 

1.9 ± 0.2 

1.2 ±0.2 

1.2 ± 0.2 

0.07 ±0.04 

0.62 ±0.09 

0.008 

2.2 ±0.9 

3.7 ± 1.3 

1.1 ±0.2 

1.6 ± 0.3 

1.2 ±0.2 

1.1 ± 0.2 

0.11 ± 0.06 

0.61 ±0.11 


0.020 

1.6 ±0.7 

2.7 ± 1.0 

1.0 ±0.2 

1.4 ±0.3 

1.0 ±0.2 

1.0 ±0.2 

0.2 ±0.09 

0.64 ±0.13 


0.040 

1.2 ±0.5 

2.0 ±0.7 

0.9 ±0.2 

1.3 ±0.3 

0.8 ±0.2 

0.8 ±0.2 

0.23 ± 0.11 

0.74 ±0.17 


0.001 

31.2 ±4.0 

44 ± 5.0 

12.7 ±0.7 

17.9 ±0.8 

12.0 ±0.6 

11.0 ±0.5 

0.003 ± 0.0008 

0.46 ±0.01 

10Mg/yr 

0.004 

25.5 ±3.4 

57.8 ±4.7 

11.3 ±0.6 

19.2 ±0.9 

11.9 ±0.6 

11.8 ±0.6 

0.06 ± 0.01 

0.63 ±0.03 

0.008 

22.0 ±2.9 

37.0 ±4.0 

10.8 ±0.7 

16.4 ±0.9 

12.0 ±0.6 

11.2 ± 0.5 

0.11 ± 0.02 

0.61 ±0.04 


0.020 

16.1 ±2.4 

26.9 ±3.3 

9.7 ±0.7 

14.0 ±0.9 

9.8 ±0.5 

9.9 ±0.5 

0.19 ±0.03 

0.64 ± 0.04 


0.040 

11.7 ± 1.7 

19.9 ±2.3 

8.6 ±0.6 

12.9 ±0.8 

8.2 ±0.5 

8.4 ±0.5 

0.22 ± 0.04 

0.73 ± 0.05 


0.001 

155 ± 9 

219 ± 11 

63 ± 1.3 

89 ± 1.6 

60 ±1.1 

53 ± 1.1 

0.003 ± 0.0004 

0.46 ± 0.006 

5OM 0 /yr 

0.004 

128 ±7 

290 ± 10 

57 ± 1.3 

96 ± 1.7 

60 ±1.1 

59 ± 1.2 

0.06 ± 0.006 

0.63 ±0.01 

0.008 

110 ±6 

185 ±8 

54 ± 1.2 

82 ± 1.6 

60 ± 1.2 

56 ± 1.1 

0.11 ±0.008 

0.61 ±0.01 


0.020 

81 ±5 

135 ±7 

49 ± 1.5 

70 ± 1.9 

49 ±1.1 

49 ± 1.0 

0.19 ± 0.01 

0.64 ± 0.02 


0.040 

59 ±4 

100 ±5 

43 ± 1.3 

65 ± 1.7 

41 ±1.3 

42 ± 1.3 

0.22 ± 0.02 

0.73 ± 0.03 


0.001 

310 ±11 

438 ± 14 

127 ± 1.5 

179 ± 1.8 

119 ±1.4 

106 ± 1.3 

0.003 ± 0.0003 

0.46 ± 0.004 

lOOMg/yr 

0.004 

252 ± 9 

574 ± 11 

113 ± 1.5 

192 ± 1.9 

120 ± 1.5 

118 ± 1.5 

0.06 ±0.002 

0.62 ± 0.008 

0.008 

219 ±7 

369 ± 10 

107 ± 1.7 

164 ± 2.2 

119 ± 1.2 

112 ± 1.0 

0.10 ±0.005 

0.60 ±0.01 


0.020 

162 ± 6 

271 ±8 

98 ± 1.6 

140 ± 2.2 

98 ± 1.7 

99 ± 1.6 

0.19 ±0.009 

0.64 ±0.01 


0.040 

119 ±5 

201 ± 7 

86 ± 1.6 

130 ±2.1 

82 ± 1.6 

83 ± 1.3 

0.22 ± 0.01 

0.73 ± 0.02 


different strong-line calibrations are presented in Table [3] 

We take T04 model as our standard metallicity classification 
and find lower metallicity galaxies dominate in number as 
they are mainly dwarf galaxies. They contribute little to 
the FUV flux and a small but significant fraction of the Ha 
flux. None of the CCSNe occur in the lowest metallicity 
galaxies. The galaxies from LMC (Z=0.008) to twice Solar 
metallicity (Z=0.04) however contribute most of the light 
and CCSNe in our sample. We use the relative fractions 
of the FUV flux at different metallicities to create our 
model population because of a better agreement between 
FUV flux and the CCSN rate. The Ha flux peaks towards 
lower metallicities and is sensitive to a shorter period of 
star formation, only the last 10 Myrs rather than the 
longer timescale for CCSN and FUV flux of 100 Myrs. 
Therefore the FUV flux better represent the averaged 
SFR during the period where SN progenitors will be born. 
We would expect the low metallicity dwarf galaxies to be 
younger with more Ha dominated star b urst than high 
metal licity galaxies. Also recent work by iBoauien et al.l 
(l2015h suggested that the amount of Ha emission that 
escapes depend on the intensity of star formation in a 
galaxy. From their result we expect more Ha emission 
to escape from dwarf galaxies but less Ha emission to 
escape from giant galaxies which have more diffuse star¬ 
forming regions. However the main reason why the Ha and 
FUV fluxes are greater is that lower metallicity stars are 
hotter and therefore emit more Ha flux relative to FUV flux. 

The variation of emitted fluxes with metallicity can 
be seen in our model populations shown in Table [4] It 


presents the predicted Ha and FUV flux at different 
metallicities from single-star and binary-star populations 
with star formations of lMg/yr, lOMgyr” 1 , SOMgyr” 1 , 
lOOMgyr” 1 along with uncertainties in the prediction. 
These fluxes are derived through our synthetic galaxy 
model which is discussed in detail in Section 3.3. The SFR 
is contributed from total galaxy’s star formation which 
vary with stochastic effect from constructing galaxy by 
individual star clusters. Therefore at given SFR the uncer¬ 
tainties of flux, CCSN rate and CCSN relative ratio, are 
highly dependent on how many star clusters are required. 
In addition, the uncertainties decrease with increasing SFR 
as stochastic effect decrease. The FUV fluxes for both single 
and binary models, decline slightly as metallicity increases 
compared to the Ha flux which changes by a factor of 2. 
Moreover, FUV flux experience similar declining as CCSN 
rate does when metallicity increase. As for the relative ratio 
of type Ibc to type II SN, it gradually increases with the 
increasing metallicity but keeps almost constant at given 
metallicity whatever the value of total SFR is. In addition, 
compared to single-star model binary-star model produced 
more Ha and FUV flux and increased the relative ratio of 
type Ibc and type II SN, but made little difference in total 
CCSN rate. This will be discussed in detail in section 3.2. 
On balance, FUV flux is our best choice for this metallicity 
classification for the 11HUGS galaxies. 

The choice of different metallicity calibrations can lead to 
large variations in the classification of 11HUGS sample as 
shown in Table [3] Compared to the T04 model, the KK04 
model is biased towards higher metallicities in all aspects: 
observed galaxy fraction, the fraction of Ha and FUV flux 
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Table 5. Ha and FUV flux, CCSN rate and ratio predictions at lMg/yr of star formation rate following different metallicity relations. 


Metallicity 

L Ha /(10 41 

erg s 1 ) 

Lfuv/(10 28 

erg s 1 Hz 1 ) 

Rccsn/( 10 2 yr) 

R'lbc/ Rll 

Model 

single 

binary 

single 

binary 

single 

binary 

single 

binary 

T04 

KK04 

PP04 03N2 

1.39 ± 0.59 
1.34 ± 0.57 

1.82 ± 0.77 

2.41 ±0.77 
2.29 ±0.79 
3.26 ± 1.06 

0.95 ±0.17 
0.96 ±0.17 
1.07 ±0.16 

1.45 ±0.30 
1.45 ±0.28 
1.62 ±0.30 

0.97 ±0.22 
0.98 ±0.21 
1.08 ±0.19 

0.98 ±0.23 
0.99 ±0.22 
1.12 ±0.18 

0.14 ±0.12 
0.15 ± 0.10 
0.14 ±0.08 

0.60 ±0.15 
0.59 ±0.15 
0.61 ±0.12 


and observed CCSN fraction. Conversely the PP04 model 
shows a preference to lower metallicities. Even though there 
are still around 70 per cent of galaxies in KK04 model laying 
below solar metallicity, over 90 percent of Ha and FUV flux 
and total CCSN are from solar and twice-solar metallicities. 
On the other hand PP04 03N2 model push more galaxies 
below solar metallicity and almost half of the total Ha 
flux, 40 per cent of FUV flux and 30 per cent of CCSN are 
given from galaxies below solar metallicity. These differ¬ 
ences in metallicity classifications are finally resulting in the 
variations of our metallicity-mixed models shown in Table [5] 

We find, compared to T04, the KK04 model favours 
higher metallicities, decreasing the Ha flux slightly and no 
effects on FUV flux for both single-star and binary-star 
models. While the PP04 model that is biased to low metal¬ 
licities, enhances the Ha emission by a third and slightly 
increases FUV emission. This result is consistent with what 
we get from a single composition as shown in Table [4] with 
more Ha and FUV fluxes from lower metallicities and larger 
variations in Ha flux than FUV flux. However different 
metallicity classifications only lead to a slight variation in 
CCSN rate and the relative ratios from different metallicity 
models even keep almost constant. This is because for all 
the metallicity classifications the CCSN rate is dominated 
by either solar (PP04 03N2 model) or twice-solar (KK04 
model) metallicity galaxies. Table []] shows a little differ¬ 
ence in the CCSN rates and the relative ratios between 
these two single metallicity models. In summary, we can 
match the observations better by changing the metallicity 
compositions of galaxies. In our following work we adopt 
T04 model as a standard to make metallicity classification 
of 11HUGS sample. These supplementary models we have 
created allow us to estimate how uncertain our predictions 
are due to assumed metallicities of the 11HUGS galaxies. 

3.2 Binary Effect on CCSN Rate and SFRIs 

We construct our galaxy model using the results of BPASS, 
a population-synthesis model based on a d etailed stellar 
evolution code calculated in lEldridge et al.l (;2 008j ). There 
are two sets of detailed evolution model, single stars and an 
extensive set of binary stars, to produce a realistic synthetic 
stellar population as observations show approximately 70 
per cent of massive sta rs have their evo lution affected by 
binary interactions (e.g. [Sana et~aL IEqH). These evolution 
models are combined with an atmosphere model to predict 
emitted spectrum of the stars. Moreover, a composite stellar 
spectrum is calculated for a range of ages. This process is 
described in detail in Eldridge & Stanway (2009, 2012). 


To calculate the CCSN rate, first we determine whether 
a star produces CCSN and which types CCSN occur at 
the end of its evolution. We use the same requirements for 
both single-star population and binary-star population. A 
star with final mass greater than 2 Mq and carbon-oxygen 
core mass over 1.38 Mq are taken to experience a CCSN. 
All the identified CCSN progenitors are then divided into 
type II SNe with hydrogen in the final model and type 
Ibc SNe without hydrogen. Therefore, the fate of a star 
at given initial mass is determined by its mass loss during 
its lifetime. Compared to single-star model, mass transfer 
in an interacting binary systems can remove the outer 
hydrogen envelope from a star’s surface and also lead 
to mass gain from the companion star. This result is a 
difference in the CCSN rate of different types as shown 
in Figure [T] The predicted rates are from the composite 
model obtained by mixing the five metallicity populations. 
For both single-star and binary-star model, type II SNe 
arise from similar stellar population in terms of stellar age 
arranging from 6 Myrs to 60 Myrs, because more massive 
stars that die earlier have strong stellar winds that will 
remove the hydrogen for both sets of models. The binary 
model however have a reduced occurrence of type II SNe 
and increased the probability of type Ibc SNe. In addition, 
a type Ibc SN progenitor can only occur after a lifetime of 
10 Myrs from an interacting binary star. When it comes to 
the total CCSN rate, binary populations make little differ¬ 
ence, but as we have demonstrated it plays an important 
role in changing the relative ratio of type Ibc to type II SNe. 

Interacting binaries also have a great effect on both 
Ha and FUV fluxes. Although the effect diverges with age. 
The Ha line from a star-forming galaxy is produced by 
hydrogen gas that have been photonionized by hot massive 
stars. This emission is c ontributed by stars with stellar 
ages less than 16 Myrs dKennicuttl [l9 98; Go garten et al.l 
l2009ll . Therefore, we derived Ha luminosity by taking 
account of all the Lyman continuum photons N(LyC) 
below the wavelength of A = 912A and foll ow t he hydrogen 
recombination relation provided by Osterbrock fe Ferland 
( 120061 ) : 

L(Ha) = 1.36 x 1CT 12 N(LyC) erg s _1 (2) 

While the FUV flux arises from the continuum spectrum 
at the wavelength between 1556A to 1576A directly from 
the stars themselves spectrum. Stars co ntribute to this 
with ages up to appr oximately 100 Myrs llKennicu ttl ll998l : 
iGogarten et ahlElOOglh 

How the Ha and FUV fluxes vary with stellar age of 
our composite metallicity population are shown in the 
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log 10 (age/yr) log 10 (age/yr) 



log 10 (age/yr) log 10 (age/yr) 


Figure 1. Distributions of CCSN rate in terms of stellar age for type Ibc SNe in the left top panel and for type II SNe in the left bottom 
panel. The relations of Ha and FUV flux across stellar age are shown in right top panel and right bottom panel separately. The red solid 
lines are from single-star model and blue dashed lines are from binary-star model. 


right two panels of Figure [T) We see that the binary-star 
population produces more flux than single-star population. 
Th ese are in ag reeme nt with the greater values predicted 
bv lZheng et al.l (l2012l '). The extra emission occurs because 
of the binary interactions creating more hot stars with 
their hydrogen envelope stripped and mass-gain in binary 
systems. The FUV fluxes follow each other with the binary 
populations being only 0.1 to 0.2 dex greater at any age. 

A greater difference occurs for the Ha flux. After 3 
Myrs w hic h is the life duration of a typical H II region 
dCrowtherjl2013li . the Ha flux from binary-star model does 
not decay as rapidly as single-star model. This again is due 
to more longer-lived stars but also more hot massive stars 
at later times to continue to contribute to the Ha flux. 
In addition, the short lifetime of the typical H II region 
implies that binary-star models have only a small effect 
on observed H II regions. However over longer timescale 
the amount of diffuse Ha emission from a galaxy might 
represent less star-formation than previously thought. 


3.3 Our Synthetic 11HUGS Model 

Using the BPASS results described above, we construct our 
synthetic 1 1HUGS ga laxie s model as follows. We follow the 
method of lEldridgg (120121 ) by creating galaxies from indi¬ 


vidual star clusters from a range of masses between 50Mq 
and 10 6 Mq. We pick the masses at random ac cording to a 
cluster mass function from lHunter et ahl (l2003l l of the form: 


dN , 2 

— oc M 2 
dM 


(3) 


Rather than uniformly distributing the age of clusters we 
set their age randomly too. For example to reproduce a 
constant star formation we allow the cluster to have any 
age from 0 to 10 Gyrs ago. With each cluster’s age and 
mass we then select and scale the appropriate BPASS 
output and add it into our total predicted spectrum and 
current CCSN rate of model galaxies. After repeated work 
as above for different metallicity galaxy populations we mix 
the model galaxies based on the metallicity classification of 
11HUGS sample and obtain a composite 11 HUGS model. 
With this method of creation we have a prediction of 
SFRIs appearance and CCSN rate from stellar populations 
spanning a SFR. from 10 _3 MQyr _1 to lOOMgyr -1 . At low 
SFRs our results do predict a range of possible appearances 
due to stochastic effects which can be seen in Table [3] 
fe.g. IPf lamm-Altenburg et. alJ 1 20091 : iFumagalli et al.l 1201 ll : 
lEldridgel 20121 : Weisz et al. 2012h . however in the amount of 
SFR of the 11HUGS sample are at a level beyond stochastic 
effects. 
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t/yr 

Figure 2. Delayed SFR as a function of time under 4 different time 
and compare with constant SFH in red solid lines. The right panel 



t/yr 


scale parameter with logarithm value of 6.3, 6.5, 7.0 ,7.5 respectively, 
is in linear time scale and the left panel is in logarithm time scale. 


A further variation to our 11 HUGS model is that we 
have also crea ted models with a delayed SFH as in 
Lee et al.l (|2010h . The SFR, 'F(Ur), in the delayed SFH is 
given by 

^(t, t) oc \e~ t/r (4) 

where t is the time since the onset of star formation and 
r is the star formation time scale parameter. In our work 
we allow r to vary from 1 Myrs to 30 Myrs. Compared to 
constant SFH, delayed SFH has its dominated star-burst 
time range around t where we can expect relatively high 
SFR. appears. This time range is directly dependent on 
the choice of r. Increasing r increases the proportion of 
older stars relative to the younger stars, and therefore 
alters the predicted results of our synthetic models. As 
shown in Figure [2] SFH with lower values of r experience 
a sharper increase at earlier times. As r increases, delayed 
SFH approach to constant SFH evolving similarly. From 
the plot in logarithm time scale, we can have a more 
clearly look at its dominated star-burst period. Here, the 
position with the largest increasing slope correspond to a 
star-burst peak which is around r. After the peak SFH 
decline gradually and a lower value of r speed up this 
decline. As s result, increasing the value of r does not 
only delay the star-burst peak time, but also extend the 
duration of the star formation. Moreover, we can associate 
our predicted CCSN rate to a stellar population with its 
dominated star-burst period known. 

Our synthetic galaxy models under a delayed SFH 
and constant SFH predicting CCSN rate and full syn¬ 
thetic spectrum can be compared directly to the observed 
11HUGS galaxies. 

4 COMPARISON WITH 11HUGS 
4.1 Constant SFH Model 

Our constant SFH model is the simplest case to consider. 
We show how the Ha? and FUV fluxes vary relative to 


Table 6. Relative ratio of type Ibc SNe to type II SNe (Ribc/Rn) 
under constant SFH. 


* 


Ribc/Rn 


Observed ratio 


0.31 ±0.18 


z 

single 

binary Mix single &; binary 

0.001 

0.002 

0.45 

0.23 

0.004 

0.07 

0.63 

0.35 

0.008 

0.10 

0.60 

0.35 

0.020 

0.20 

0.65 

0.42 

0.040 

0.23 

0.72 

0.48 

Mixed Z 

0.22 

0.70 

0.46 


the CCSN rate in Figure [3] The values expected for the 
11HUGS galaxies are shown by the boxes indicating the 
range of values allowed for the observed fluxes and CCSN 
rates. We find both single-star and binary-star models 
perform reasonably well with predictions within a factor 
of 2 to the observed values. Also the agreement is closer 
for the FUV flux rather than the Ha flux. However the 
models tend to overestimate the amount of Ha and FUV 
fluxes required to reproduce the observed CCSN rate, 
showing there is a p ossible excess of CCSNe as found by 
several studies fe.g. I And o et al.1 I 200!tI: iKistler et all l201ll : 

iMattila et al.l 12012 : 1 11 o r i u ch i et al.1 2013fl . We also find that 

the standard single-star population achieves the closest fit 
to the observed parameters. Our binary-star models pro¬ 
duce much greater Ha flux than the single-star population 
which causes tension with the observations. However there 
are limitations to our models which currently do not yet 
account for the effects of dust on the Ha and FUV emission. 

On the other hand there is another test we can apply 
as shown in Table [6] The ratio of the rate of type Ibc 
SNe to the rate of type II SNe is 0.31 ± 0.18. The large 
error is due to the small number of events. This ratio is 
strongly dependent on age, interacting binary fraction and 
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Cumulative CC SNe Rate per year Cumulative CC SNe Rate per year 


Figure 3. Constant star foramtion model: Hct and FUV flux as a function of cummulative CCSN rate. The red solid lines are from 
single-star model and blue dashed lines are from binary-star model. The colored boxes show the observed flux and CCSN rate of Sample 
A, B, C as shown in Table [T] Since the FUV flux of Sample A is arichived from our estimation, we mark it in a dashed red box. 


metallicity. The effect of the latter is shown in Table [6] 
The predicted CCSN rate ratios for our single-star and 
binary-star models are either side of the observed value. 
The single-star prediction is lower than the observed value 
but within the range allowed by errors. By contrast, the 
ratio from binary stars is higher than the observed value. 
However a mix of single-star and binary-star population 
also gives a good fit to the observed value. 

4.2 Delayed SFH Model 

With constant SFH we predict too much Ha and FUV 
flux relative to the CCSN rate. We therefore decided 
to investigate alternative SFHs in order to reduce the 
contribution of young stars. A delayed SFH is applied as 
described in Section 3. The outcomes we achieve under the 
delayed star birth assumption are shown in Figure [4] Here 
the shaded regions show the range of Ha and FUV fluxes 
from our synthetic galaxies in which the predicted CCSN 
rate matches the observed rate of the 11HUGS galaxies.We 
see that for matching the Ha flux and the CCSN rate at 
the same time the preferred values for log(r) are between 

6.3 to 6.7, or times of 2 to 5 Myrs. Although the binary-star 
models are only just below the observed values and so 
longer values of r would not be unreasonable. While to 
match the FUV flux and CCSN rate a much wider range of 
possible values for log(r) are possible between 6.3 to 7.4, or 
times of 2 to 25 Myrs. This implies that FUV flux follows 
on a much longer timescale than the Ha flux as expected. 

Figure [5] illustrates how the relative ratio of the two 
CCSN type rates, Ribc/Rn, varies with r. This plot is 
similar to Figure 0] showing that as t increases the binary- 
star models are favoured relative to the single-star models 
where the ratio drops for the single stars. The key difference 
between the single-star and binary-star models is that the 
single-star population can only reproduce the observation 
for r in a short period of approximately 2 Myrs. While the 
binary population, or a mix of single stars and binaries 
could reproduce the 11HUGS observations over a much 
longer time scale. However the delayed SFH models are in 


effect just a starburst model where the age of the burst is 
linked to its duration. So the best fit we can rule out is that 
the SFH that best matches the 11HUGS galaxies is one 
that is not dominated by a recent young burst but older 
more mature bursts. 

In summary, our delayed SFH models show a prefer¬ 
ence for 3 Myrs in matching with observed values. This 
fact could be due to a sample selection effect but this is 
remote as 11HUGS are a distance limited sample not Ha 
selected. It is more likely that we are over predicting the 
Ha emission over the first 3 Myrs of the star formation. 


5 DISCUSSION & CONCLUSIONS 

We have created a synthetic model of the 11HUGS 
galaxies to investigate the effect of interacting binary 
stars. We find binaries do not vary the overall CCSN 
rate to any great degree but do boost the number of type 
Ibc SNe and the Ha and FUV SFRIs linked to hottest stars. 

In trying to match the observations of the 11HUGS 
galaxies we must discuss the caveats within our population 
models. We have not yet accounted for the absorption 
of the ionizing photons by dust grains and therefore 
these photons do not form Ha flux (e.g. iRelano et al.l 
l2012h . This effect is not yet included in our models. In 
addition, the 11HUGS samples are only nearly complete. 
As in Sample A there is more star formation occurring 
than in Sample B, the expected Sample A FUV flux 
might be greater than we predict here. This result would 
then be in closer agreement with the binary-star popu¬ 
lation under constant SFH model. Therefore we should 
not be too surprised that the match is not perfect but 
there is more work to refine the synthetic galaxy population. 

The most surprising fact however is that our constant 
SFH models produce too ma ny CCSNe in t h e nea rest 11 
Mpc of space as argued in iHoriuchi et al.1 (I2013T ). This 
implies that all supernovae predicted by the models have 
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29.6 
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Figure 4. Delayed SFH model: Ha and UV flux as a function 
of timescale parameter r. The shaded regions show the range of 
Ha and FUV flux when the predicted CCSN rate matches the 
observed rate of the 11HUGS galaxies in Table [T] The pairs of 
black straight lines show the observed value range of Ha and UV 
flux in Sample AB and C as shown in Table [T] Because the FUV 
flux of Sample A is arichived from our estimation, we mark it by 
two dashed black lines. 


to be observed. If some massive stars form black holes at 
core collapse a nd are unobservable event as suggested by 
iHoriuchi et alj d2014h then this excess of stellar deaths 
compared to stellar population models becomes more 
severe. However estimating whether a specific stellar model 
would produce an observable SN or not requires detailed 
modelling of the core-collapse and the subsequent explosion. 

On the other hand, the delayed SFH model argues 
that the predicted CCSN rate from Ha and FUV SFRIs is 
dependent on the assumption of the SFH. One interpreta¬ 
tion of our preferred delayed SFH model is that it implies 



log 10 (x/yr) 


Figure 5. The relative ratio of type Ibc SNe and type II SNe 
(Rlbc/Rll) as a function of timescale paparamter r. The red dot¬ 
ted line is from sinlge-star model and blue dotted line is from 
binary-star model. A mixed population of single-star and binary- 
star model is shown in the yellow dotted line. The green horizontal 
lines show the observed CCSN ratio range (0.31 ±0.18) as shown 
in Table [6] 


we over estimate the flux from young stars and therefore 
get an excess of CCSNe. A similar result is also provided 
by ionizing photon leakage and dust attenuation on Ha 
and FUV flux. Therefore to truly determine if the SFRIs 
and CCSN rate match, the ionized photon leakage and 
dust attenuation must be considered. However this is not 
straightforward. 

We have also investigated how sensitive our results 
are to the metallicity calibration we use for modelling 
the metallicities of the 11HUGS galaxies. We find that 
having the more metal rich biased population leads to small 
changes in the Ha and FUV fluxes. With a more metal 
poor distribution we would find even more Ha and FUV 
fluxes and would result in greater disagreement between the 
observations and our models. Therefore the mismatch be¬ 
tween our constant star-formation models and observations 
is difficult to explain by an incorrect metallicity distribution. 

Our main conclusion is that single-star models matches 
the observations best in terms of the prediction of total 
CCSN rate, when the effects of ionizing photon leakage 
and dust attenuation are not considered and assuming the 
local star formation has followed a constant SFH over the 
last few 100 Myrs. Binary-star models however prolong the 
stellar age of stars to go CCSN and therefore provide a 
better prediction of CCSN rate under various SFH spanning 
star-burst duration between 3 and 25 Myrs ago. In addition, 
single-star models under constant SFH predicted a lower 
ratio of type Ibc relative to type II SN but still within 
the observed uncertainty. By contrast, binary-star models 
overestimate the ratio. However, a mix of singl e-star an d 
binary-star model also gives a good fit as in lEldridee et al .1 
(|2008i . 120131 ). The limitation in making a firm determination 
is due to the large uncertainty caused by limited numbers of 
CCSN discovered. With more CCSN discovered within 11 
Mpc this uncertainty will decrease. If the number of CCSN 
is doubled we estimate that it may become possible to rule 
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out the single-star only population. One final argument 
against the single-star model is that under the delayed SFH 
the single-star model only match the observations for a 
burst of star formation only 2 Myrs old. However mixing in 
binaries to the single-star population provides a better fit 
with weaker constraints on the age of the population. 

Finally we note that our results are robust if the CCSN 
sample is complete. However we already obtain too many 
CCSNe compared to its associated Ha and FUV. If there 
are missing SNe or other massive stellar deaths in our 
sample that are not included in our CCSN sample then the 
disagreement between models and observations becomes 
worse with too many stars dying compared to the number 
we observe forming. 
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